Static dielectric response and Born effective charge of BN nanotubes from ab initio 

finite electric field calculations 
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Ab initio investigations of the full static dielectric response and Born effective charge of BN nan- 
otubes (BN-NTs) have been performed for the first time using finite electric field method. It is 
found that the ionic contribution to the static dielectric response of BN-NTs is substantial and also 
that a pronounced chirality-dependent oscillation is superimposed on the otherwise linear relation 
between the longitudinal electric polarizability and the tube diameter (D), as for a thin dielectric 
cylinderical shell. In contrast, the transverse dielectric response of the BN-NTs resemble the be- 
havior of a thin (non-ideal) conducting cylindrical shell of a diameter of D + 4A, with a screening 
factor of 2 for the inner electric field. The medium principal component Zy of the Born effective 
charge corresponding to the transverse atomic displacement tangential to the BN-NT surface, has 
a pronounced D-dependence (but independent of chirality) , while the large longitudinal component 
Z* exhibits a clear chirality dependence (but nearly D- independent), suggesting a powerful way to 
characterize the diameter and chirality of a BN-NT. 



PACS numbers: 73.63.Fg, 77.22.-d, 78.67.Ch 
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I. INTRODUCTION 

Since the discovery of carbon nanotubes (CNTs) in 
1991 1 , carbon and other nanotubes have attracted con- 
siderable interest worldwide because of their unusual 
properties and great potentials for technological appli- 
cations. For example, semiconducting CNTs could func- 
tion as nanoscale field-effect transistors 2 . It was also pre- 
dicted that nanotori formed from metallic CNTs may ex- 
hibit giant paramagnetic moments, 3 Though CNTs con- 
tinue to attract great interest, other nanotubes such as 
BN nanotubes (BN-NTs)^ are interesting in their own 
right and may offer different possibilities for technologi- 
cal applications that CNTs cannot provide. In particular, 
as far as the optical and opto-electronic applications of 
nanotubes are concerned, BN-NTs could be superior to 
CNTs because BN-NTs are uniformly insulating, inde- 
pendent of their chirality^ Furthermore, recent experi- 
ments indicate that BN-NTs exhibit stronger resistance 
to oxidation at high temperatures than CNTs. 6 

Therefore, the properties of both single-walled and 
multiwalled BN-NTs have been intensively studied ex- 
perimentally and theoretically in recent years. In this pa- 
per, we focus on the dielectric response, including static 
dielectric constant, electric polarizability and Born ef- 
fective charge, of BN-NTs. A detailed knowledge of di- 
electric response is needed, e.g., to characterize optical 
excitations, screening at contacts, and plasmons in nan- 
otube arrays. Variations in size and chirality may or may 
not affect dielectric properties, which in turn may be ex- 
ploited to align nanotubes during, e.g., plasma-enhanced 
chemical vaposition synthesis 7,8 and to separate different 
tubes in solutions^. Synthesis and separation of specific 



nanotubes remain an important challenge. 

Previously, the static dielectric constant and electric 
polarizability of nanotubes were typically obtained as the 
zero frequency limit of the optical conductivity (see, e.g., 
Refs. [lOliill and references therein), which is calculated 
using either semiempiral tight-binding method or ab ini- 
tio local density approximation (LDA) or generalized gra- 
dient approximation (GGA) within independent particle 
linear response theory (IPLRT). In this way, only the 
electronic contribution to the static dielectric constant 
(eoo) is calculated] 10 ' 1:L i 12 ] 13 Furthermore, local field ef- 
fects, which can be substantial for the electric field ap- 
plied perpendicular to the tube axis?^, are neglected] 10 ! 11 
Here, we apply the latest finite-field electric-enthalpy 
theor y 14 : 15 to study the static dielectric properties of BN- 
NTs. This theory allows us not only to obtain the full 
dielectric constants including both electronic and ionic 
contributions but also to include the local field effects. 
Furthermore, also using this theory, we perform the first 
systematic study of the Born effective charge of the BN- 
NTs. A knowledge of the dependence of Born effective 
charge on the diameter and chirality of the nanotubes not 
only is important to understand their dynamical proper- 
ties but also would allow us to characterize the nanotubes 
using, e.g., Raman and infrared spectroscopies. 



II. THEORY AND COMPUTATIONAL DETAILS 



We use the latest finite field electric-enthalpy 
theoryi^ii 5 . to calculate field-induced electric polariza- 



tion P and atomic forces F. For 
E, the dielectric constant is e a p — 



a small finite field 

S a p + -lirXa/3 where 
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Xa/3 = P a /Ep. If the atoms are kept fixed, this yields 
the electronic contribution (so-called clamped ion di- 
electric constant) and if both the electrons and atoms are 
allowed to relax in response to the field, the full static di- 
electric constant e s t a tic is obtained. The electric polariz- 
ability a is given by e a p = S a p + 47ra a /3 / fl where is the 
unit cell volume. The Born effective charge on the jth 
atom is eZ* Q/3 = Fjp/E a . Traditionally, the Born effec- 
tive charge is obtained via eZ* a „ — Q(dP a / 'dujp)\^ = o 
where Uj is the displacement of the jth atom due to 
the optical phonon mode, and is usually calculated by 
the finite difference approach. Therefore, the finite field 
calculations have an advantage of being simple and less 
CPU-time consuming especially for complex systems like 
BN-NTs. 

We consider a number of representative BN-NTs with 
a range of diameters from all three types of BN-NTs, 
namely, the zigzag [(n,0), n = 5, 6, 8, 12, 16], armchair 
[(n,n), n = 3, 4, 6, 8, 12], and chiral (4,2), (5,2) BN- 
NTs. For comparison, we also study the dielectric prop- 
erties of bulk BN in the cubic, wurtzite and hexagonal 
structures (denoted as c-BN, w-BN and h-BN, respec- 
tively) as well as a single honeycomb BN sheet. Our ab 
initio calculations were performed using accurate projec- 
tor augmented-wave (PAW) method^, as implemented 
in the QMAS (Quantum MAterials Simulator) package^. 
A supercell geometry was adopted so that the nanotubes 
are aligned in a square array with the closest interwall 
distance between adjacent nanotubes being 9 A. The iso- 
lated BN sheet is simulated by a slab-supercell approach 
with an inter-sheet distance of 15 A. A plane- wave cut-off 
of 35 Ry was used throughout. For BN-NTs, a 2 x 2 x n 
mesh in the first Brillouin zone is used throughout. For 
the achiral nanotubes, n — 8 and for the chiral nan- 
otubes, n = 6. The ideal nanotubes were first con- 
structed by rolling- up a honeycomb BN sheet. Their 
atomic positions and lattice constants were then fully 
relaxed until all the atomic forces and the stress were 
less than 4.0xlCT 5 Ry/Bohr and l.OxlCT 6 Ry/Bohr 3 , 
respectively. The same force and stress convergence cri- 
teria were used in the finite-field calculations. The small 
finite-field used is 0.001 atomic unrfe!£. The calculated 
field-induced polarization and Born effective charge for 
the field perpendicular to the tube axis were checked 
with the usual saw-tooth potential method which can- 
not be applied for the field parallel to the tube axis, and 
the results from both approaches agree with each other 
very well (up to the third digit). Note that the unit cell 
volume f2 is not well defined for nanotubes and the single 
BN sheet. Therefore, like the previous calculations^!!, 
we used an effective unit cell volume of the nanotubes 
and the BN sheet rather than the volume of the super- 
cells which is arbitrary. The effective unit cell volume 
for a nanotube and the BN sheet is given, respectively, 
by £1 = TiDdT and by f2 = Ad where d is the thickness 
of the BN sheet which is set to the interlayer distance 
of /i-BN. D and T are the diameter and length of trans- 
lational vector of the nanotube!!, respectively. A is the 



TABLE I: Lattice constants (a,c), static dielectric constant 
(e), and Born effective charge (Z*) of bulk BN in cubic, 
wurtzite and hexagonal structures (denoted as c-BN, w-BN 
and h-BN, respectively) as well as the isolated BN sheet. The 
available experiment values are listed in brackets. 
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area of the surface unit cell of the BN sheet. Further- 
more, in the supercell approach to a single wire (an iso- 
lated slab), the periodic dipole images can give rise to 
a substantial depolarization field which not only affect 
the transverse component of the field-induced properties 
but also decrease rather slowly as the supercell dimen- 
sion increases!^. Therefore, the transverse dielectric con- 
stant, electric polarizability and Born effective charges 
presented below, have been corrected for this artificial de- 
polarization field via one-dimensional Clausius-Mossotti 
relation!^ for the single BN sheet and two-dimensional 
Clausius-Mossotti relation! 3 -!^ for the BN-NTs. We have 
checked numerically that the transverse component of 
these physical properties are independent of the super- 
cell dimension used (within 1 %) for the single BN sheet 
and also the (5,0) BN-NT. 



III. RESULTS AND DISCUSSION 

A. Bulk BN structures and single BN sheet 

The results for the bulk BN structures which are sum- 
marized in Table I. All the theoretical lattice constants 
of all bulk BN structures agree rather well (within 1.5 
%) with the experimental values (see Table I), except 
c of h-BN which is 3.3 % smaller than the experimetal 
value. Interestingly, the static dielectric constants and 
Born effective charges of c-BN and w-BN are very simi- 
lar, and the anisotropy for w-BN is rather small (Table 
I). This is due to the fact that although the layer stacking 
along the c direction of the wurtzite structure and the one 
along [111] direction of zincblcnde structure are different, 
the local tetrahedral atomic configuration is common to 
both structures. The calculated dielectric properties for 
c-BN and w-BN agree rather well (within 4.0 %) with 
the available experimental values (see Table I). In con- 
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FIG. 1: (Color online) (a) Atomic arrangement of a single BN 
sheet and (b) Schematic illustlation of the mode of deforma- 
tion of BN sheet under the in-plane electric field. Note that 
the magnitude of deformation is exaggerated. 



trast to w-BN, the anisotropy in the dielectric properties 
of h-BN is rather pronounced. For example, the Born 
effective charge for the field along the c-axis (Z c ) is only 
about I /3 of that for the field perpendicular to the c-axis 
(Z a (b)) due to weak covalency along the c-axis. Here a 
denotes one of the two in-plane primitive vectors (i.e., 
the vector connecting the two nearest neighoring B(N) 
atoms) and b is the vector that is perpendicular to a. 
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FIG. 2: (Color online) (a) Longitudinal (a zz ) and (b) trans- 
verse (ctxx) polarizabilities of BN-NTs as a function of tube 
diameter (D). The dashed lines in (a) are the straight lines 
fitted to the calculated polarizabilities (symbols), and in (b) 
are the function of b(D + S) 2 fitted to the calculated values. 



It is interesting to compare the last two columns (h-BN 
and BN sheet) in Table I. As h-BN is a one-dimensional 
stacking of BN sheets with rather large separation be- 
tween neighboring sheets, it is natural to expect that 
both systems may have similar in-plane properties, as 
is actually seen in Table I except the static dielectric 
constant. While there is no in-plane anisotropy for the 
clamped ion dielectric constant and Born effective charge 
for both h-BN and BN sheet, only the static dielectric 
constant of BN-NT shows strong in-plane anisotropy. 
The reason for this is that lattice strain is induced in 
a BN sheet under the in-plane electric field. The mode 
of strain is schematically shown in Fig. 1, where the 
maginitude of strain is much exaggerated. The strain 
mode is clearly anisotropic with respect to the direction 
of electric field and induces the anisotropy in the in-plane 
static dielectric constant for BN sheet. However, such a 
lattice strain cannot be induced in h-BN because of the 
alternating array of B and N along the c-axis. The B and 
N arrangement of a given sheet is reversed in the neigh- 
boring sheets so that the strain induced by electric field 
in the neighboring BN sheets cancels out. 



B. BN-NT 

1. Longitudinal dielectric constant and electric 
polarizability 

The calculated static dielectric polarizability a zz and 
dielectric constant e s t a tic of the BN-NTs are plotted as 
a function of diameter (D) in Figs. 2 and 3, respec- 
tively. Fig. 2(a) shows that on average, both the full 
and clamped ion (electronic contribution only) polariz- 
abilities a zz per unit length for the electric field parallel 
to the tube axis (E \\ z) is roughly proportional to the 
tube diameter, i.e., a = aD, resembling the behavior of 
a thin dielectric cylinderical shell. This behavior is in 
strong contrast with the latest finding of a zz ~ D 2 + C 
for the CNTs from the density-functional perturbation 
theory (DFPT) calculations^ and also with the recent 
finding of from the IPLRT calculations^. Because the 
number of atoms on a BN-NT per unit length is pro- 
portional to D, this result indicates that every atom on 
the BN-NTs has nearly the same static polarizability, as 
would be the case for insulators. The slope of a zz is a z 



= 4.57 A and a z ,o = 3.08 A for the full and clamped ion 
polarizabilities, respectively. Note that the slope a z $ is 
nearly identical to the previous IPLRT calculations (3.04 
Aii, 3.13±2) indicating that local field effects for E \\ z 
are indeed negligible. However, when the ionic contribu- 
tion is included, as in the present calculations, the electric 
polarizability a zz is significantly increased, by ~ 50 %, 
demonstrating the importance of the ionic contribution. 

There is clear chirality dependence in a zz and e zz , as 
shown in Fig. 2(a) and Fig. 3. We first discuss the 
clamped ion dielectric constants which are denoted by 
inverse triangles in Fig. 3 to simplify the problem. For a 
given diameter of BN-NT, zigzag BN-NTs have larger di- 
electric constant than armchair ones and the chiral ones 
are intermediate. The chirality dependence is strong for 
smaller diameters and tends to vanish as the diameter in- 
creases. The longitudinal dielectric constant of clamped 
ion case clearly converges to 4.7 which is the clamped 
ion dielectric constant of a single BN sheet, which is ex- 
pected to be the case. The strong chirality dependence of 
clamped-ion longitudinal dielectric constant comes from 
the chirality dependence of the band gap. It was re- 
ported in Ref. |ll| that the band gap of the zigzag BN- 
NTs would be reduced substantially as the tube diameter 
decreases, while that of the armchair BN-NTs depends on 
the diameter only slightly. It is known that the systems 
with a smaller band gap would have a larger dielectric 
constant. For static longitudinal dielectric constant, the 
present calculation takes account of only the relaxation 
of atomic positions with a fixed unict cell. Due to the 
difficulty in practical calculations, the shear deformation 
of the unit cell as observed in a single BN sheet (Fig. 1) 
is not taken into accout. Therefore the behavior of the 
static longitudinal dielectric constant shown by triangles 
in Fig. 2 is basically the same as those of clamped ion 
case. The large diameter limit of static dielectric con- 
stant is about 6.4 which is close to the value for a single 
BN sheet obtained within the same constraint of atomic 
relaxation. 



2. Transverse dielectric constant and electric polarizability 

We begin our discussion with rather simple features. 
First, the transverse quantities (in Fig. 2(b) and Fig. 
3) do not show appreciable chirality dependence for the 
whole range of diameter. Second, the ionic contribution 
is rather small, being only ~ 15% of the electronic con- 
tribution in the polarizability [Fig. 2(b)]. Third, the 
large diameter limit of transverse dielectric constant may 
be given by the circular average of the response of each 
part of the BN-NT surface, which may be approximated 



by 0.5(e 



aa,bb 
static 



.) with quantities for a single BN 



sheet. As we do not take account of the shear strain 
for BN-NT, e"^ tfc is about 6.4 and the converged value 
of the squares (in Fig. 3) in the large diameter limit 
may be about 4.0. The most surprising fact in the trans- 
verse quantities is that the transverse polarizability a xx 
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FIG. 3: (Color online) Full static dielectric constant of BN- 
NTs as a function of tube diameter. Each dashed line is the 
curve fitted to the corresponding calculated dataset. Longi- 
tudinal clamped-ion dielectric constants of BN-NTs are also 
plotted for comparison. 



in Fig. 2(b) exhibits a rather different behavio r 11 ' 13 . In 
fact, a xx (D) can be nicely fitted to a xx = b(D + S) 2 , as 
can be seen from Fig. 2(b). This indicates that as far as 
the transverse dielectric properties are concerned, a BN- 
NT behaves like a thin (nonideal) conducting cylinderical 
shell^ (or a system of electrons moving freely on a cylin- 
drical shelliS) with an effective diameter of D+5. The pa- 
rameters b and 6 obtained by the fitting are, respectively, 
0.0629 and 4.13 for the full polarizability, and 0.0523 and 
4.40 for the clamped ion polarizability. 6 is slightly larger 
than the effective thickness (d ~ 3.3A) of the BN sheet 
mentioned above. The screening factor of the inner elec- 
tric field (E in ) by a BN-NT is E out /E in = 1/(1 - 86) w 2 
(E out is the applied field outside the BN-NT) when the 
full response of the BN-NT is considered. Note that for 
a long classical dielectric cylindrical shell, the inner elec- 
tric field Ej n is also uniform and along E ollt . However, 
E in = 4E out e/ [{e+l) 2 -(e-l) 2 (D + d) 2 /(D-d) 2 } would 
be enhanced rather than screened. For example, for the 
(8,8) BN-NT, the screening factor is 0.4 if e x s f atic (2.58) 
[for the (8,8) BN-NT] is used or 0.8 if ef taUc (1.65) [for 
the BN sheet] is used. To see whether this phenomeno- 
logical interpretation in term of a thin conducting cylin- 
drical shell is an accidental coincidence or not, we plot 
the calculated electric potential on a plane perpendicular 
to the tube axis of the (8,8) BN-NT in Fig. 4(a) and also 
calculate the total electric field along the applied electric 
field line going through the center of the BN-NT [Fig. 
4(b)]. Fig. 4 clearly shows that the electric field inside 
the BN-NT is indeed uniform and considerably reduced. 
The screening factor from Fig. 4(b) is 2.08 but not in- 
finite as would be expected from a perfect conducting 
cylindrical shell. We also plot the electric field for some 
other BN-NTs and find the same situation. Note that 
both the metallic and semiconducting CNTs exhibit the 
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FIG. 4: (Color online) Upper panel: Potential contour plot on 
a xy-plane perpendicular to the tube axis. Lower panel: The 
electric field along the line parallel to the :r-axis and going 
through the center of the BN-NT in the xy-plane shown in 
the upper panel. 



N 



~i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — r 

(5,(^(6,0) « j ( i2 0) (16,0) " 



(4,: 



- V -V— v — -V- 



(3^3) (4,4) (6,6) (8,8) 



■ D z i 

A Z* 
V Z* 



A-A- 



A-A 



.,-A 



---A" 



fc-B -E3-- — B - £3— □ — -B - 



10 



12 



14 



Diameter (A) 



FIG. 5: (Color online) Three principal Born effective charges 
of BN-NTs as a function of tube diameter. Each dashed line 
is the curve fitted to the corresponding calculated dataset. 



same behavior with a larger screening factor of 4.4J£ 
Therefore, the transverse polarizability from the present 
finite field calculations reveals a qualitatively different 
behavior, when compared with previous IPLRT calcula- 
tions ; 11 ! 13 demonstrating the importance of taking the 
local field effects into account in this case. 



3. Born effective charges 

To obtain the Born effective charges, we perform for 
each BN-NT three self-consistent finite-field calculations 
for the field applied along three Cartesian coordinate 
axes, respectively. A Born effective charge tensor for each 
atom is then constructed from the field-induced atomic 
forces. These Born effective charge tensors are then diag- 
onalized to give three principal components (Z*, Z* Z*). 



Interestingly, it turns out that Z* , Z* are for the trans- 
verse fields perpendicular to and tangential to the local 
tube surface of the ion concerned, respectively, while Z* 
corresponds to the field along the tube axis. The calcu- 
lated principal Born effective charges on all the BN-NTs 
studied here are displayed as a function of the tube di- 
ameter in Fig. 5. 

Figure 5 shows that the medium component (Z*) of 
the BN-NTs exhibits a pronounced monotonic depen- 
dence on the tube diameter (but independent of chirality) 
while the small component (Z*) is almost a constant. In 
contrast, the large component (Z*) shows a clear chi- 
rality dependent oscillation superimposed on the mean 
value of 2.70 e (the in-plane components of the Born ef- 
fective charge of the single BN sheet) (see Fig. 5 and 
Table I). These interesting diameter-dependence of Z* 
and chirality-dependence of Z* suggest that one could 
determine both the diameter and chirality of a BN-NT 
by measuring its Born effective charges using, e.g., Ra- 
man and infrared spectroscopies. Characterization of the 
chirality and diemeter of a BN-NT is still a great chal- 
lenge. 

Interestingly, as the tube diameter becomes large, Z* 
and Z* tend to approach to the Z* and Z*^ components 
of the Born effective charge of the single BN sheet for the 
field perpendicular and parallel to the sheet plane (Table 
I), respectively. This is perhaps not surprising because 
the observed monotonic diameter dependence of the Z* 
and Z* merely reflect the degree of the local curvature of 
the BN-NTs with a finite diameter. 



IV. CONCLUDING REMARKS 

To summarize, we have carried out a systematic ab ini- 
tio study of the static dielectric properties of the BN-NTs 
within density functional theory with the LDA plus the 
latest finite-field electric-enthalpy theory. Specifically, 
the full static dielectric constant, electric polarizability 
and Born effective charges of a number of BN-NTs as 
well as the single graphitic BN sheet and bulk BN struc- 
tures have been calculated. The calculated lattice con- 
stants, dielectric constants and Born effective charges for 
the bulk BN structures are in good agreement with avail- 
able experiments. We find that the ionic contribution to 
the static dielectric response of BN nanotubes (BN-NTs) 
is substantial and that a clear chirality-dependent oscila- 
tion is imposed on a linear relation between the longitu- 
dinal electric polarizability and the tube diameter, as for 
a thin dielectric cylinderical shell. We also find that the 
calculated features of the transverse dielectric response of 
the BN-NTs resemble that of a thin conducting cylindri- 
cal shell, with a screening factor of 2 for the inner electric 
field. Our calculations also show that for each BN-NT, 
the medium principal component of Born effective charge 
tensor has a pronounced dependence on diameter (but in- 
dependent on chirality) and the large principal componct 
exhibits a chirality dependence (but nearly independent 
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on diameter), suggesting that a BN-NT could be char- 
acterized by measuring its Born effective charges using, 
e.g., Raman and infrared spectroscopies. 
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